CLAIMS 



1. Method for estimating the seismic illumination 
fold \{x,p) in the migrated 3D domain at at least one 
image point x / for at least one dip of vector p 

wherein . the illumination fold I(x,p ; s,?) for each (source- s r 
receiver F ) pair in the seismic survey is estimated, by 
applying the following steps: 

- determination of the reflection travel time 
t r ( x r (p) ; J , F ) from the source s to the specular 
reflection point x r on the plane reflector ' passing 

through the image point x and perpendicular to the dip 
vector p and then returning to the reflector F ; 

starting from the diffraction travel time t d (x;s 3 r ) 
from the source s to the said image point 3c and then 
returning to the reflector F ; 

- incrementing the said illumination fold I(x,p;s,?) 
related to the said (source s, receiver F) pair as a 
function of the difference between the diffraction travel 
time t d {x;s,r) and the reflection travei time 

tr(i r (?);i;r). 

2. Method according to claim 1, comprising the step 
of summating each of the said illumination folds Iu.p;!,?) 
related to a (source s, receiver F) pair so as to 
determine the total illumination fold I(x 9 p) = yKx.^iJ) . 

s,p 

3. Method according to one of the preceding claims, 
wherein, during the incrementing step, the illumination 
fold I(x,p,s,r) is incremented using an increment function 
i(t d ,t r ; s,r ) according to I(x,p) = I(x,p) + i(t d/ t r ; s,r), the 
said increment function taking account of the difference 
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between the diffraction travel time t d (x; s, r) and the 
reflection travel time t r (x r (j5); s, r). 

4. Method according to claim 3, wherein the 
increment function i is a function of the seismic wavelet 
s(t). 

5. Method according to claim 4, wherein the 
increment function i is expressed as a function of the 
derivative of the seismic wavelet s(t) according to: 
i(t d , t r ; J, r) = s(t d (x; s , r) - t r (x t {p);s,r) 

6. Method according to claim 4, wherein the 
increment function i is expressed as a function of the 
derivative 5 (t) of the seismic wavelet s (t) with respect 
to time according to: 

i(t d , t r ; s , r) = t d (x; s , r) - t r ( x r ( p ) ; s , r ) 

1. Method according to any one of claims 3 to 6, in 
which an a priori correction w(x,s,r) of the illumination 
fold is taken into account by migration, comprising the 
step of incrementing the illumination fold I{x,p;s,r) 
related to a (source s, receiver r) pair by 
i (t d , t r ; s \ r )..w(x; s , r ) . 

8 . Method according to any one of the preceding 
claims, wherein the determination step includes the 
second order Taylor series development of the diffraction 
travel time t d {x;s,r) around the image point 3c: 

t«(X;S , r )=t d (X; :«,?) + ( V x t d (X;S , r )) T . (X r -X) + - ( x r - 
X) T .A s , x t d (X; S ,r ).(X r -X) 

9. Method according to claim 8, wherein the specular 

reflection point x ' ( p ) is determined along the length of 
the said reflector such that the diffraction travel time 

at the said specular reflection point x »- (P ) is 
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stationary, according to the equation: 

p t A( V x t d (x ; s , r )+A x , K t d (* ; s,r ) . (X r (P)-X)) = 0 . 

10. Method according to any one of claims 8 or 9, 
wherein the specular reflection point x r and the 

5 reflection travel time t t {x s [p);s,r) (are determined 
according to the following expressions: 

x r (p) = x-M.F" 1 . b 

t r (x r {p ) ; s , r ) = t d ( x ; s , r ) . b r . F" 1 . 6 

10 

where : 

- M is a (3x2) matrix described by two vectors 
extending along the length of the reflector, and 
therefore perpendicular to the dip vector p ; 

15 - b is a (2 x 1) vector of first order derivatives 

of the diffraction travel time along the reflection 

plane: b = M T . ( V x t d ) ; 

" - F is a (2 x 2) matrix of second order derivatives 
of the diffraction travel time along the reflection 
20 plane: F = M T . (A x#3C t d ) .M. 

11. Method according to any one of the preceding 
claims, wherein the determination step uses isochronic 
migration maps t d {x;s,r ) specified for each (source s , 
receiver F ) pair involved in the migration at each image 

25 point x in the migrated 3D. domain. 

12. Method according to any one of the preceding 
• claims, wherein the seismic illumination fold I{x 9 p) in 

the migrated 3D domain is estimated during the Kirchhoff 
summation migration of seismic data recorded during the 
30 3D seismic prospecting. 
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13- Method for correction of seismic data amplitudes 
recorded during 3D seismic prospecting in order to 
compensate for the effect of non-uniform illumination of 
sub-soil reflectors, comprising the steps of: 

- estimating the illumination fold I(x,p) using the 
method according to any one of claims 1 to 12, 

- using the inverse r x {x 9 p) of the said ratio as a 
weighting factor . to be applied to each of the said 
seismic data amplitudes. 

14. Method for selection of an acquisition geometry 
among a plurality of acquisition geometries as a function 
of the target of 3D seismic prospecting, comprising the 
steps of: 

- determining the illumination fold I{x,p) by the 
method according to any one of claims 1 to 12, for each 
of the acquisition geometries considered, 

- selecting the acquisition geometry providing the 
optimum illumination fold as a function of the target. 
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